•v 



«1» > 




NAVAL POSTGRADUATE SCHOOL 
Monterey, California 




THESIS 

MAGNETIC FIELDS GENERATED BY 
INTERNAL OCEAN SEAWATER MOTION 

by 

Charles A. Davis 
December, 1991 

Thesis Advisor: Dr. Donald L. Walters 

Approved for public release; distribution is unlimited. 



T260069 



UNCLASSIFIED 

SECUPi"^^ CLASS^FtCAT Qr, QP ^his PAr g 



REPORT DOCUMENTATION PAGE 



form Approved 
OM8NO 0704 0188 



la REPORT SECURITY ClASSiFiCATlON 

UNCLASSIFIED 



ib RESTRICTIVE fVAPK ,\GS 



2a security CLASSiFiCAT OfvJ AUTHORITY 



2b DECLASSIFICATION UOWNORAD NG SCHE DUlE 



3 D'STRiB'J^'ON Av/AlAB _TY OF RtP-'^"^ 

Approved for public release; distribution 
is unlimited. 



4 PERFORMING ORGANIZATION REPORT NUMBER(S) 



5 MONiTORiNG ORGANiZATiON REPORT N^rV'BER;S; 



HP 






6a NAME OF PERFORMING ORGANIZATION 



6b OFFICE SYMBOL 
(If applicable) 



7a NAME OF MON TORiNG ORGAN ZAT 0^ 



6c ADDRESS {City, State, and ZIP Code) 



7b ADDRESS ^C/fy State and ZIP Code) 



8a NAME OF FUNDING SPONSORING 
ORGANIZATION 



8b OFF’CE symbol 
(If applicable) 



9 PROCUREMENT INSTRUMENT iD£NTiF;C A'^ON N UMBER 



8c ADDRESS (Gty, State, and ZIP Code) 



’'■) SO..PCE OF F lND'NG N 'MBtPS 



PROGRAM 


PROJECT 


TASI' 


vVOP-< jNlT 


Element no 


NO 


NO 


ACCESSION NO 



11 TITLE (Include Security Classification) 

Magnetic Fields Generated by Internal Ocean Seawater Motion 



12 PERSONAu AUThOR(S) 
Davis, Charles A. 



13a TYPE OF REPORT 

Master’s Thesis 



13b T.ME COVERED 
FROM TO 



date of report {Year, Month Day) 

December, 1991 



IS PAGE COuNT 

80 



16 supplementary notation The views expressed in this thesis are those of the author and do 
not reflect the official policy or position of the Department of Defense or the 
U. S. Government. 



17 COSAN CODES 


FIElD 


GROUP 


SUBGROUP 















18 Subject terms {Continue on rt^erse if necessary and identify by block number) 

Induced magnetic fields; Magnetic anomaly detection; 
Internal Waves; Boundary layer turbulence 



19 abstract {Continue on reverse if necessary and identity by block numoer) 

This thesis models induced magnetic fields from the motion of seawater in the earth’s 
magnetic field analytically and compares the results to arctic on-the-ice magnetic 
fluctuation measurements. The oceans have various types of internal motions, such as 
internal waves and turbulence. This motion of seawater, which is a conductor, in the 
earth’s magnetic field induces a current density. This current density, in turn, 
induces its own magnetic field. This thesis models internal waves and upper layer 
ocean turbulence analytically. The corresponding induced magnetic fields are 
calculated using a static form of Maxwell’s equations and parameters for the Arctic 
are inserted. Comparisons are made with measurements from the Arctic Internal Wave 
Experiment (1985). The predicted fields from internal waves have magnitudes that are 
measurable and of the same order of magnitude as ionospher ically generated fields. The 
predicted fields from turbulence are several orders of magnitude smaller than 
iotiobpherlcaily; geLteiated fields. besides giving information about internal motions in 



20 DISTRIBUTION AVA LABNiTY Of ABSTRACT 

[jj UNCuA SS'F'ED ' iNLiMi^ED □ SAME AS RP^ □ DT|C USERS 



21 ABSTRACT SECuRiTy C.ASS F-CATiON 

TINrT.AF;qTFTFr) 



22a NAME OF RESpoNLiB^E ^NDlvlDuAu 

D. L. Walters 



^b telephone (/nc/uUe Area Code) 

408-646-2267 



22c QFFiCE SYMBOu 

PH/We 



DD Form 1473, JUN 86 



Previous editions are obsolete 

S/N 0102-LF-014-6b()3 



SEULRMr ClAb^.F_CA/ 



AGE 



UNCLASSIFIED 



1 



U NCLASSIFIED 

SECURITY CLASSIFICATION OF THIS PAGE 



Block 19 (cont.) 

the ocean, the seawater induced fields are a noise source in magnetic anomaly 
detection . 



DO Form 1473. JUN 86iReverse) 



SECURI'’’'y classification of this PAGf 



11 



UNCLASSIFIED 



Approved for public release; distribution is unlimited. 



Magnetic Fields Generated by Internal Ocean Seawater Motion 

by 

Charles A. Davis 
Lieutenant, United States Navy 
B.S., Duke University 

Submitted in partial fulfillment 
of the requirements for the degree of 

MASTER OF SCIENCE IN PHYSICS 

from the 



NAVAL POSTGRADUATE SCHOOL 



December 1991 



r - ■ - - 

Karlheinz E. Woehlcr, Chairman 
Department of Physics 



iii 




///^/^ 

lJ 



ABSTRACT 



This thesis models induced magnetic fields from the motion of seawater in the earth’s 



magnetic field analytically and compares the results to arctic on-the-ice magnetic fluctuation 
measurements. The oceans have various types of internal motions, such as internal waves and 
turbulence. This motion of seawater, which is a conductor, in the earth’s magnetic field induces 
a current density. This current density, in turn, induces its own magnetic field. This thesis 
models internal waves and upper layer ocean turbulence analytically. The corresponding 
induced magnetic fields are calculated using a static form of Maxwell’s equations and parameters 
for the Arctic are inserted. Comparisons are made with measurements from the Arctic Internal 
Wave Experiment (1985). The predicted fields from internal waves have magnitudes that are 
measurable and of the same order of magnitude as ionospherically generated fields. The 
predicted fields from turbulence are several orders of magnitude smaller than ionospherically 
generated fields. Besides giving information about internal motions in the ocean, the seawater 
induced fields are a noise source in magnetic anomaly detection. 
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I. INTRODUCTION 

Seawater is a conductor situated in the earth's magnetic 
field. Motion of this conductor from surface waves, internal 
waves, turbulence, etc. produces a motional electromotive 
force. The resulting current density induces its own magnetic 
field. Understanding the production and propagation of such 
magnetic signatures provides information about the underlying 
seawater motion. Also, these induced magnetic fields are a 
noise source for other magnetic measurements, such as magnetic 
anomaly detection. 

This thesis constructed analytical models of these induced 
fields and compared the predicted values with experimental 
measurements. One model considered internal wave sources and 
another ocean upper boundary layer turbulence. The internal 
wave model had versions for the Arctic and lower latitude 
oceans . 

The internal wave model started with the mathematical 
derivation of a governing equation for seawater velocity from 
fluid mechanics. The model assumed no rotation, zero 
viscosity, incompressibility, and slow variation of mean 
density with depth. Separation of variables provided a linear 
differential equation for the amplitude of the vertical 
component of velocity. The corresponding solution and the 
separation of variables equations together provided an 
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analytical model of seawater velocity, v. The velocity 
frequency dependence was then tailored to fit observed 
velocity spectra for either the Arctic Ocean or lower latitude 
oceans . 

Next, the induced magnetic field was calculated. The 
motion of the conducting seawater with velocity v in the 
earth's magnetic field B (assumed constant) induces an 
electromotive force and causes a current density J. This 
current density J induces its own magnetic field B'. Because 
the frequencies involved are on the order of 10'^ Hz or less, 
displacement currents were neglected and the inductive field 
B' calculated with the Biot-Savart Law. The total induced 
field B' at a field point was calculated assuming an internal 
wave of infinite y-extent with wavenumber in the x-z plane. 
Integration over all x and y, and over z within a depth range 
from Dq to D yielded B*. The wavenumber dependence was 
removed by two different methods. In one, B* was converted to 
a power spectrum, transformed from a one-dimensional to a two- 
dimensional spectrum, and integrated over wavenumber from a 
minimum, to infinity to give B'^((i)). In the other, 

application of a boundary condition to the vertical component 
of velocity produced a dispersion relation between w and 
Choosing the lowest mode of oscillation, using the dispersion 
relation, and converting to a power spectrum yielded B'^(w) . 

Data from the Arctic Internal Wave Experiment (AIWEX) in 
1985 provided a means to compare the model B'^ power spectral 
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density as a function of frequency with experimental results. 
The predicted magnitudes for the induced fields were in 

the same range as extrapolated magnetometer measurements. The 
large spatial extent of internal waves produced signals of the 
same order of magnitude as ionospherically generated signals. 
The first k-dependence method yielded frequency responses that 
do not follow the 1/f^ dependence of the data. The second 
method does yielded the 1/f or 1/f^ frequency dependence 
expected. The success of the second method implied that the 
modal structure of the internal wave field must always be 
utilized in calculating the induced magnetic fields. 

The turbulence induced magnetic field model had a similar 
formulation. It assumed the region of turbulent seawater 
resides in the upper boundary layer of the ocean between 
depths and D, and to be three dimensionally isotropic 
within this region. The seawater velocity v was analytically 
modeled as a linear superposition of plane waves propagating 
in three dimensions with various wavenumbers and frequencies. 
This motion in the earth's magnetic field B resulted in the 
current density J, which was utilized in the Biot-Savart Law 
to calculate the induced magnetic field B' at a field point, 
again neglecting displacement current. The total field was 
found by integrating over all x and y, and over z between 
and D. Keeping the principal term of the B' expression, the 
induced field components were integrated over wavenumber from 
a minimum, to infinity. Then the horizontal field 
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magnitude, B'„, was converted to a power spectral density as 
a function of frequency. 

This power spectrum was compared to the AIWEX magnetometer 
data. The frequency dependence of the model agreed with the 
observed 1/f^ dependence of the data. However, the predicted 
magnitudes were several orders of magnitude below the data. 
The AIWEX signal in the lO’^-lO’^ Hz range appeared to arise 
from fields generated in the ionosphere. 
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II 



BACKGROUND 



A. OCEAN ENVIRONMENT 

The conductivity of seawater principally derives from 
dissolved salts. In general, conductivity decreases with 
depth. An average value for ocean conductivity is 4 (S/m) . 

Seawater density varies in the ocean. Temperature and 
salinity are the main governing factors. Density as a 
function of depth generally displays three zones (Gross, 
1971) . In the surface zone, from the surface to 50-100 m, the 
density remains approximately constant with depth due to 
mixing of the upper layer of ocean by wave action. The next 
zone, from 50-100 m to around 500-1000 m, forms the pycnocline 
zone. Density increases rapidly with depth. Generally, the 
density follows a monotonic increase because the temperature 
declines with depth. The bottom zone, below the pycnocline, 
is called the deep zone. Here, density slowly increases with 
depth due to decreasing temperature. Fig. 1 (Gross, 1971) 
depicts these zones and their general temperature, salinity, 
and density variations. 
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Figure 1; Typical Ocean Thermocline, Halocline, and 

Pycnocline (Gross, 1972) 

The earth's permanent magnetic field permeates the oceans. 
Since the field is approximately that of a magnetic dipole, 
the direction and magnitude of the field vary across the 
oceans. It goes from basically horizontal with respect to the 
surface of the earth and magnitude of about 30,000 gamma 
(nanotesla) near the equator, to basically vertical with 
magnitude of about 60,000 gamma at the north and south 
magnetic poles. Seawater motions of the ocean within this 
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field induce the electromotive force and the resulting current 
density J. 

B. INTERNAL WAVES 

A change in density with depth ( dp/dz < 0 ) allows 

buoyancy oscillations within the ocean known as internal 
waves. The momentum and continuity equations for fluids are 
used to derive an equation for the vertical velocity of these 
internal waves. 

Following Gill (1982), the fluid momentum equation is 

-^+20xu=-p’’^Vp-g’ + vV^u, (1) 



where u = fluid velocity = ui + vj + wk, n = system rotation 
angular velocity, and v = kinematic viscosity. Assume 
isentropic motion, i.e. no viscous effects (v=0) and no 
rotation (n=0) . The mass equation of continuity is 

p-i^ + V-u = 0. 

Dt 



Assume density is a function of potential temperature, 0, and 
salinity, s, independent of pressure, p. Using the chain rule 
and taking 0 and s constant with respect to time over a 
differential element 

Dp _ dp DQ ^ dp Ds _ Q (3) 
Dt de Dt ds Dt 
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The continuity equation then becomes the incompressible fluid 
condition 

V-u = 0. ( 4 ) 



Assume small perturbations in pressure, p', and density, p'. 



P = Po + P'> P = Po + P'' 



^Po 

dz 



= -99 o' 



( 5 ) 



The fluid momentum equation becomes 



P = -Vp' - p'p^. 



( 6 ) 



Applying the convective derivative with average velocity v, 
the time derivative of density becomes 



-Pp - , 

Dt dt 



( vV) V 



dp' 

dt 



+ w- 



^Po 

dz 



0 . 



( 7 ) 



The fluid velocity components u and v can be eliminated by 
first taking the partial derivative with respect to time of 
the continuity equation 



_3 Vu . . 

dt dtdx dtdy dtdz 



( 8 ) 



and then using Equation 6 in component form to replace time 
derivatives of u,v with spatial derivatives of p'. The 
result is 
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d^w _ d^p' ^ d^p' 
^ ° dzdt dx^ dy^ 



( 9 ) 



Combining Equation [ 7 ] and the time derivative of the z- 
component of Equation [6] yields 



ac2 






j__av 

p Q dzdt ' 



( 10 ) 



where N is the Brunt-Vaisala frequency, given by 



7/2 



g ^Po 

Po cfz ■ 



( 11 ) 



Taking the horizontal Laplacian of Equation [ 10 ] 



82 , 82 



82 



8^2 dx~ dy‘ 



) w + 7/2 ( 



82 , 82 



8^2 



)ur = -^( 
p o 



82 . 82 



8x2 0y2 



82pU 2) 

dzdt ' 



Substituting Equation [ 9 ] in the right hand side and 
rearranging yields an equation for w, the z-component of 
velocity, 



82 

8t2 



[ 



8x2 Qy2 dz ° dz 



772 ( 



82 , 82 

8x2 0y2 



) w = 0 . (13) 



Z, 



Assume that the vertical velocity, w, varies with depth, 
much more rapidly than the density, Then 



1 8 , ^ \ 82 w 

p o 8z ° dz az2 



(14) 
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which is known as the Boussinesq approximation. 
Equation [13] then simplifies to 



at2 



JL+_§L 

dx^ dy^ 



+ ] w + N'^ { 

dz^ 



dx^ dy^ 



) w 



0 . 



( 15 ) 



This equation describes the behavior of the vertical velocity 
component of an internal wave disturbance. 

Separation of variables can simplify Equation [15] 

w{x,y,z,t) = Xix) Yiy) Z{z) T(t) . ( 16 ) 



The X, Y, and T portions have the form 



d^x 

dx^ 



-H 



0 . 



( 17 ) 



with the following solutions 



Z = 






y = y^e 



ikyY, 



T = T^e 



lU>t 



( 18 ) 



The z part becomes 



d^Z 

dz^ 



+ 



klil- 



N{z)^ . 
0>2 



Z = 0. 



( 19 ) 



This equation possesses the Sturm-Liouville form. Thus, its 
solutions are orthogonal eigenfunctions Z(z) with real 
eigenvalues. The vertical velocity solutions then take the 
form 
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( 20 ) 



Since they are orthogonal, any internal disturbance can be 
expanded as a sum of these solutions. 

For Z(z), try solutions of the form 



The term with N(z)^ makes this, in general, a nonlinear 
differential equation. 

To simplify this equation, approximate p(z) as a piecewise 
continuous function over adjacent regions such that dp/dz is 
a constant in each region. Then N(z) will be a constant in 
each region. Equation [22] can be solved for each region, and 
the solutions joined with boundary conditions specifying 
continuity of velocity. Alternately, approximate p(z) by a 
superposition of linear slope segments each with dp/dz equal 
to a constant. Then N(z) is a constant for each segment and 
Equation [22] is again readily solved. These approximations 
are seldom applicable to the real ocean (Gill, 1982) for 
accurate velocity profiles. However, the resulting solutions 



Z{z) = V {z) ^ 



( 21 ) 



Substituting this into Equation [19] gives for w(z) 
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do have the general shape and the oscillating character of the 
profiles resulting from actual p(z) distributions. 

With N(z) equal to a constant, Equation [22] becomes a 
second order, constant coefficient, linear differential 
equation 



dz^ 



+ a 



dz 



+ Pv = 0, 



(23) 



where a - 2iJc^, p - [ ( ) {k^ + k^) - Jc|] . 



to' 



The solutions of the characteristic equation are 



X. . . -Ik, ± 1 

± 2 z 



\ 



( ikl . k}) ( 25 ) 

(0^ 



Define 5 as 



6 



N 



, - Oi 

0)2 



2 

-) 



ikl + ky) . 



(26) 



In a stratified ocean, the maximum oscillation frequency 
equals the buoyancy frequency N. Thus, o < N , and S is real. 
For w < N, the general solution of equation [23] is 

i?(z) = e'-‘'*^''[q + Cj e-"’»^] . (27) 
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Initial and boundary conditions determine the constants Cl and 



C2. 

The vertical velocity component is then 



w(x,y, z, t) 



v(z) 



v(z) 



(28) 



The incompressible continuity condition. Equation [6], 
determines the horizontal components, u and v. 



u(x,y,z, t) = v] { 



) e 



i (Jt 



0) t) 



(29) 



v(x,y,z,c) = [i-^ - Jc v] ( — (30) 

kl . kl 

These analytical expressions for the seawater velocity now 
facilitate calculating the induced magnetic fields. 

C. ARCTIC UNDER ICE SEAWATER TURBULENCE 

In the arctic, turbulent seawater flow exists from below 
the ice (effectively zero depth) down to some depth, D, at 
which the motion transitions to internal waves. 40 m is a 
typical value for D (private communication with T. Stanton) . 

Utilize superposition of three-dimensional plane waves of 
varying wavenumber and angular frequency, o, to model this 
turbulent region 
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v= v,(;c„o>)e'**"jf ] e'^“^ (31) 



Assume three-dimensional isotropy of turbulence to simplify 
the analysis 

v^(k^,(xi) = Vyiky,oi) = v^{k^,<ji) * vik.u)) . (32) 



Empirically model the velocity amplitude function v(k,o)) as 



v(ic,(o) = 



v„ 



^Jkl + k^ ^o>o + (0^ 



(33) 



The power spectrum for this function is then 



v^{k,(j>) = 



{kl + k'^) (o>o + co^) 



(34) 



The constant factors and in the denominator keep the 

power finite as k and u> approach zero. For large k, this 
exhibits the same k'^ dependence as for internal waves 
(Garret, Munk, 1972) . For large o, this exhibits the 

dependence characteristic of under ice internal wave fields 
(ONR Report, 1991) . The velocity expression is 

v(k,(o,r,t) = 

sjkl + + ( 0 ^ 
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Ill 



MAGNETIC SIGNATURES 



A. PROPAGATION OF INDUCED ELECTROMAGNETIC FIELDS 

This model neglected attenuation due to the propagation of 
the electromagnetic waves through seawater, a conductor, 
because of the low frequencies involved and the corresponding 
long wavelengths. The electromagnetic waves are characterized 
by their wavelength, X. For a "good conductor", like 
seawater, and a "good insulator", like air, the wavelengths 
are 



X 



conductor 



2n 



X 



insulator 



C 

f 



(36) 



where 



|i o = cond= 4 — , c=3 .0x10^— . (37) 

m s 

Table I lists wavelength values for typical frequencies of 
interest in internal ocean seawater motions. 
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Table I: Wavelength Values 



Freq (Hz) 


^seawater 


^alr . 


10° 


1570 


3x10° 


10*'' 


5030 


3xl0’ 


10*2 


15,700 


3xl0'° 


10'^ 


50,300 


3xl0" 


10*^ 


157,000 


3x10^2 



The induced magnetic fields will "propagate” with 
characteristic wavelength ^sea„ater through the ocean and then 
with the much, much longer wavelength above the surface. 

Characteristic ocean buoyancy frequencies are around 10'^ - 
10'^ Hz. Since internal wave frequencies are always less than 
the buoyancy frequency, and the depth of the Arctic Ocean is 
around 3000 m (Dietrich, et al., 1980), the corresponding 
wavelengths of the induced electromagnetic radiation are much 
larger than the depth. Within the ocean and in the air above 
the ocean surface, conduction current, J, and not displacement 
current, edE/dt, creates the induced magnetic field. The 
field point is in the "near field" and propagation effects are 
negligible. Relevant turbulence frequencies range around 10’^- 
10'^ Hz. We confined our interest to turbulence in the upper 
boundary layer of the Arctic Ocean. Since a nominal depth for 
this layer is 40 m, the electromagnetic radiation wavelengths 
involved are again much larger than the depth and fields above 
the surface are "near field" with negligible propagation loss. 
Thus, we neglected attenuation due to propagation in the 
seawater conducting medium. For these same reasons, the 
sea/ice/air interfaces do not attenuate the induced fields. 
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The earth's magnetic field inside the ocean will possess 
the same value as in the air above the ocean. This follows 
from the assumption that the earth's field is a static (DC) 
field, and that the magnetic susceptibilities of the air and 
ocean are the same (i.e. both are non-magnetic) . The 
equivalence derives from the resulting boundary conditions 
between the two media 

^ sea, tangential ~ ^aiz, tangential’ ^sea, normal ~ ^air, normal' 

B. INTERNAL WAVE SIGNATURE 

1. Geometry 

The method of Crews and Futterman (1962) allowed us to 
calculate magnetic signatures. The current density followed 
from the seawater velocity in the earth's magnetic field. 
This current density was then utilized in the Biot-Savart 
method of calculating magnetic induction at a field point. 

Fig. 2 shows the geometry of the earth's magnetic field 
within the ocean (northern hemisphere) . To simplify the 
calculations, choose the wave vector in the x-z plane. Then 
the y-components of the wave vector and the velocity are zero. 
Assume the wave extends infinitely in the ±y directions. The 
seawater in the volume element dxdydz at the point (x,y,z) 
with velocity v(x,y,z,t) forms a source element for an induced 
magnetic field at the field point located at (0,y,-h). Fig. 
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3 shows this geometry, "z" represents the depth below the sea 



surface. "h” represents the height of the field point above 
the sea surface. 




▼ 

z 



X 



Figure 2: Geometry of the Earth's Magnetic Field 
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R = X + (h+z) 



Figure 3: Geometry for Calculation of Internal Wave Induced 

Magnetic Field 
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2. Induced Magnetic Field 

Seawater behaves as a relatively good conductor, with 
conductivity of approximately 4 S/m. The motion of this 
conductor in the earth's magnetic field creates an 
electromotive force. This in turn leads to a current density 
given by 



Since velocity is confined to the x-z plane, the current 
density J is in the y-direction. 

This current density creates a secondary induced magnetic 
field. The Biot-Savart law gives the differential induced 
magnetic field due to a current element di 



The prime indicates the secondary induced field. Integration 
of Equation [40] gives the components of the total induced 
field. The assumption of infinite extent of the wave in the 
±y directions (i.e. symmetry about y=0) results in the y- 
component summing to zero. The x and z components are 



*7=0 { V X B) = a { /. 



( 39 ) 




di = J dA. 



( 40 ) 




( 41 ) 




( 42 ) 
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Substituting in the expressions for di, v, B, and R and 
integrating [see Appendix A] gives 









( k,-i6)z 

— )cosycos4> ♦ 



-^ein4> ( 









) ] 



D 



( 43 ) 









[ (- 



-k*id 



<-Jr,-id)x • 

tt — )cosycos 4> ♦ -^ein4>( 

k*io ic_ 



-k^^ib 
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( 44 ) 



Here, S reduces to 



b 



0)^ 



( 45 ) 



These equations are fairly complicated. However, they 
simplify somewhat when evaluated specifically for conditions 
in the Arctic 



B = B ]c; 4) = 90"^; COS(}) = 0; sin4> = 1; D. = 0 



( 46 ) 



The induced magnetic field components become 









[ 1 .( 47 ) 



-k^^ib 



k^-^ib 



k^iN^-2oi^) 



B[=iJi|2££e-‘«*e-^-, 



^ . 1 20 )^ ] .( 48 ) 
u’ -A:,*i6 k^*i& k^{N^-2<j>^) 



Simplifying, 



iB^=B^=(n„BoOe~*»^e-^-“^( - -- ^ )4-I K- e'*'"(-;r^cos6a - SsinfiB) ].(49) 

kl*b^ ' ' 



Inserting the expression for S, these become 
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[1 - ©"**^(003 






ic^- 



o>^ 



sin^ 






Note that the x and z components are exactly the same except 
for the phase factor of i 



Bz = i 



wi th i 




(51) 



Thus, the X and z components are 90“ out of phase, and the 
induced magnetic field vector rotates in the x-z plane. It 
has magnitude 

|B1 = |b'| = |B^| . (52) 

Important characteristics of B' are (1) exponential decrease 
in magnitude with height, h, above the internal wave field, 
and (2) increase in magnitude with the vertical height D of 
the internal wave field. 

These expressions contain angular frequency u, wavenumber 
k^, and time t. Comparison with experimental data often 
requires the total field power spectrum as a function of 
frequency. Thus, all of the wavenumber and time dependencies 
must be removed. 

The frequency and wavenumber dependencies of the velocity 
were empirically modeled based on actual measurements of 
velocity power spectra. The amplitude of the internal wave 
velocity, C, which so far has only been a multiplicative 
component in the equations, became a function of frequency and 
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wavenumber. This, in turn, added additional frequency and 
wavenumber dependence to the induced magnetic field B' . 

The Arctic Internal Wave Experiment (AIWEX) , conducted in 
1985 in the Canada Basin of the Arctic Ocean (depth > 3500 m) , 
measured horizontal velocity with an electromagnetic current 
meter at a depth of 100 m. Fig. 4 shows one such spectrum. 
Sampling rate limited the highest frequencies measured to 
around 6 cycles per hour (1.7x10*^ Hz), the local buoyancy 
frequency. The spectra indicate a 1/f frequency dependence 
over the range 10'^ - 10 cycles per hour (2.8x10'^ - 2.8x10'^ 
Hz) (Levine, et al.,1987). This frequency dependence was 
empirically modeled by 

2 2 

{"a" for AIWEX) . (53) 

y/o)o + 0)2 2n^fl + 

For o >> 0 )^, it behaves as 1/u. For o < its value stays 

approximately constant and limits the total energy in the 
spectrum as o approaches zero. 

The Coordinated Eastern Arctic Experiment (CEAREX) in 1989 
measured velocity pov;er spectra at depths of 50, 100,150,200, 
and 250 m below the surface of the ice. The 50 m measurements 
were made with an acoustic doppler profiler and the other 
measurements were made with electromagnetic current meters. 
Appendix C contains reproductions of the resulting velocity 
power spectra. All plots showed an approximately 1/f^ 
frequency dependence in the range 0.1-10 cycles per hour 
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Frequency (cph) 



Figure 4: AIWEX Horizontal Velocity Spectra (solid lines) 
Garrett-Munk Model Spectrum (dashed lines) for comparison 
(Levine, et al., 1987) 



24 



(2.78x10'^ - 2.78x10'^ Hz). Often this dependence flattened 
out for frequencies below about 4x10'^ cycles per hour 
(1.11x10'^ Hz). (Czipott and Podney, 1989). 

This frequency dependence was empirically modeled by 

= — = — , {"C" for CEAREX) 

+ f2) 

For (i) >> (i)^, it behaves as l/w^, and for o < it remains 

approximately constant. Choosing corresponding to 4x10'^ 
cycles per hour gives 

0 )^ = 2nf^ = 2n(l.llxI0-’ Hz) = 1 .OxlO''' rad/ sec. (55) 

Measurements of internal wave spectra in other oceans of 
the world at lower latitudes have indicated a characteristic 
1/f^ frequency dependence. Garrett and Munk (1971) formulated 
a popular empirical model based on these observations with 
such a frequency dependence. However, the AIWEX spectra 
revealed less total energy in the arctic internal wave field 
by factors of 15-30 than in fields described by the Garrett- 
Munk model (Levine, et al.,1987). The CEAREX data indicated 
a wavefield less energetic than lower latitude fields, but 
still significantly more energetic than the AIWEX wave field. 
The CEAREX data, due to location, were subject to large 
diurnal tidal motions which pumped energy into the internal 
wave field. In addition, under ice topography (ice keels) and 
local sea floor geography ( i.e. Yermak plateau) affected the 
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internal wave field (Stanton) . In contrast, the AIWEX data 
were taken in the Canada Basin of the Arctic Ocean (Padman, et 
al.,1990) and experienced far more quiescent conditions. 

Arctic internal wave fields possess different 
characteristics from lower latitude fields. Some suggested 
causes in the Arctic are less forcing due to wind, a damping 
effect caused by the ice cover, lack of surface waves, and 
weak large scale circulation (Levine, et al.,1987) . Choice of 
frequency dependence in the internal wave field model 
critically depends on the location that is to be modeled. An 
arctic model will use a 1/f dependence for the velocity power 
spectrum, whereas a lower latitude model will use 1/f^. The 
CEAREX location represents a special situation where large 
tidal forcing causes the internal wave field to possess 
characteristics similar to lower latitude fields. 

Garrett and Munk (1971) show a collection of data on 
displacement power spectra versus wavenumber obtained from 
towed measurements. The data follow a 1/k^ pattern over a 
full four decades. The largest wavenumber corresponds to a 
wavelength of 10,000 m. Assuming a simple plane wave 
description, the velocity spectrum will possess the same 1/k^ 
dependence 

(56} 

displacement = , velocity = ^ ^^_C_QUkx-ut) ^ 
sfp- yfp- 

A simple empirical model for the k dependence could resemble 
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the omega dependence model with a power limiting factor. 
However, to facilitate later mathematical computations, use 
instead 



|2 ^ 




< k < oo. (57) 



The total energy in the spectrum is now limited by the cutoff 
at instead of having the spectrum level out at low k. 

Assume this model applies to both arctic and lower latitude 
internal wave fields. 

Combining both forms, a model for the velocity spectra for 
the Arctic (i.e. AIWEX) is 

' ^in< k < CO. (58) 

^0) o+(0^ 



The model for lower latitude and CEAREX velocity spectra 
becomes 



kcl" = 



Cc 






(0)o + k^ (fl + f^) k^ 



Ka^< k < 00.(59) 



Comparison with the AIWEX and CEAREX data determine the 
constants of proportionality, C^ and C^., At 1 cph, the AIWEX 
spectrum is approximately 0.3 (cm/s)Vcph (0.11 (m/s)^/Hz) , and 
all the CEAREX spectra have approximately the value 0.5 
(cm/s)Vcph (0.18 (m/s)VHz). Approximating the AIWEX model 
as proportional to 1/f at this frequency 
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( 60 ) 



|2 ^ 



z^r J.dk= 

2n 2n 



Ci 1 



f Kin' 



Evaluating at 1 cph 



\Vj^{l cph) P = 0.3 = 






cph 2n (1 cph) 



( 61 ) 



Then 



C^ = 2n(0.3) icm/s)^ = 1.9x10-* {-^) . ( 62 ) 



Approximating the CEAREX model as proportional to 1/f^ at this 
frequency 



I'^c 



cl 






r J_ 



dk = 



Cc 



4n^f^ Kin ' 



( 63 ) 



Evaluating at 1 cph 



|vc(l cph) p = 0.5 
' ^ ' cph 



A TT 2 / *1 2 



( 64 ) 



Then 

C| = 4nM0.5) {cm/s)^-cph Km = 5. 5x1 O'® Km - ( 65 ) 

The formulas for the induced magnetic field involve the 
velocity amplitude and not the velocity spectral density. The 
conversion to a velocity amplitude was made using the 
definition of the power spectral density. The conversion for 
this case was simple because of the simple time dependence. 
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Given a function f(t), the power spectral density is the 
Fourier transform of the autocorrelation function of f(t). 
Let F(f) be the Fourier transform of f(t). The correlation 
theorem says (Brigham, 1988) 

r- i66) 

autocorrelation of fit) =j f ix) f it+x) dx ^ FT pair - F{f)F"ir). 

The power spectral density is thus F(f)F*(f). Assuming the 
velocity has the following time dependence 

vit) • fit) = C (67) 

then the Fourier transform of the velocity is 

Vif) = J” C = C . (68) 



The power spectral density becomes 



|U(f)P = Vif) V* if) = b^if-f^) = C^. 



( 69 ) 



Thus, for this simple time dependence and real coefficient C, 
the power spectral density is obtained simply by squaring the 
coefficient. 

Applying this result to the above model, the velocity 
becomes 



v.ik,u>) = 



(w 



>)l/4 y. 



v^ik,<ji) = 









( 70 ) 



+ 
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This simple model in which the frequency dependence factor and 
the wavenumber dependence factors are multiplicative factors 
makes the following calculations much easier. However, as 
Garrett and Munk (1971) point out, actual observations 
indicate that a given wavenumber's contribution depends on the 
frequency. Following Garrett and Munk, the function n(<a) was 
included here as a reminder of this additional frequency 
dependence 



v^(k. Cl)) 



^ (w) 

(ci) o + (i)^) k 



v^ik.u)) 



Cc (ci>) 
^Ci)o k 



Kin< 



Equation [71] replaces the constant "C” in the induced 
magnetic field expression. Equation [50]. 

Examine two different methods to remove the k-dependence . 
For the first method, assume internal waves exist with a 
continuum of frequencies from zero to the buoyancy frequency, 
N, and assume they possess a continuum of wavenumbers from a 
minimum, k^^jn/ to infinity. The rationale is that actual 
internal wave fields result from a complex superposition of 
waves generated by different sources propagating in a non- 
uniform ocean. The assumption of a continuum may circumvent 
a detailed description of this wave field. The second method 
utilizes the eigenfunction mode structure of the wave field 
required by Equations [17] and [19]. Application of a 
boundary condition on the vertical velocity at depth D yields 
a dispersion relation between o and k^ for each mode. These 
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dispersion relations remove the dependence in B' . The 
power spectral density B'^ is calculated for the lowest order 
mode n=l. The lowest order modes will contain a majority of 
the energy in an internal wave field. 

Try the continuum approach first. To simplify 
computations, assume that the depth D >> Then the last 
two terms in the induced magnetic field expression, which are 
weighted by exp(-k^D) , can be neglected with respect to 1, 
simplifying the conversion to power spectral densities. The 
induced magnetic field has the simple e’“* time dependence 
referred to above. Squaring the amplitude gives the 
corresponding power spectral density 







1 (cij) (0^ + 

kx (<0 o +( 0 ^) 



(73) 



The formulas above are "one-dimensional" formulas in that 
they were derived assuming only one-dimensional wavevectors, 
k^. The more general wavefield will be two-dimensional with 
wavevectors of magnitude 

k = sjkl^kl . (74) 

The one-dimensional expressions for the induced magnetic field 
power spectrum must be converted to a two-dimensional form. 
Yaglom (1987) shows a way to accomplish this. Given a one- 
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dimensional power spectral density function, f,(k,). Then the 
corresponding two-dimensional power spectral density function 
is 



f^(k) 



2^ p d fi(iCi) dk^ 

n J, dk, yifrfff ■ 



( 75 ) 



The one dimensional induced magnetic field k dependence was 
1/k^^. Appendix D contains a summary of the details in 
computing a two dimensional k dependence. Then 

( 76 ) 



Now, perform the integration over k < k < <») 




3 1 
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dk = - 



3 1 

16 Jc4 
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k*- 



( 77 ) 



The induced magnetic field two-dimensional power spectral 
density becomes 



| b '|2 = 



16 kl 



(fa)) (N~-<x>^) (0" 






2 , ^ 2 
o+CO 



( 78 ) 



l^cl^ = 



(fa)) 



(AT^-fa)2) 0)^ 



16 kL 



N*{oi 0+^^) 



( 79 ) 



Now apply the second method using the modal structure. 
Combining Equations [26], [27] and [28] gives the vertical 
component of velocity 

w= (qe"'»^ + C 2 e"^*^) (for ky=0) . ( 80 ) 
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Let C, = -Cj = C/2i. Then 

Csin(6z) (81) 

As a boundary condition, require that the vertical velocity 
equal zero at the lower boundary of the internal wave field at 
depth D. Then 

Csin(dP) - 0, (82) 



sin (bD) = 0 , 



bD = 












= rm , n=l,2 . . , 



(83) 



= 



mi 
D 



0)" 






= 



klD'^N^ 



n^Ti^ + kln^ 



(84) 



Replace in the formula for B' in Equation [50] with this 
dispersion relation. Taking h=0 and mode number n=l gives 

-n/ 3 T 

B' = H„oBCe-^'“'^ — ^(Ar2-o)2) [ 1 + e V ] . (85) 

° tiN^ 

Converting to a power spectral density 

2 ^ ^ 

= [n jjOB] ^ — (AT^-O)^)^ [ 1 + e V ]2 (86) 

ti^N^ 

Model the velocity spectra with only a frequency dependence 
since a dispersion relation exists relating k^^ to w, 
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(87) 
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Replace in the expression, Equation [86], with Equation 
[87] 






Cl D'^ 
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(88) 



l^cp = 



0 ) 



cl g2 

2+0)2 Tc2i\^4 



(^2-o>2)2 [ 



1 + 



]2. (89) 



C. UNDER ICE TURBULENCE MAGNETIC SIGNATURE 

1. Geometry 

The same geometry as for the internal wave calculation was 
used, except that now the y-coordinate was important (Fig. 5) . 
The distance R was given by 

|R| = \jx'^ + y2 + (h+z)^ , with R = x£ + yj + z^. (8®) 

For calculational purposes, turbulence was assumed to be 
confined to the region D^<z<D. 

2. Turbulence magnetic signature 

The expressions for the earth's magnetic field and the 
induced current density are 
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Figure 5: Geometry for Calculation of Turbulence Induced 

Magnetic Field 
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B = B ic, and J = a {v X B) . 



( 91 ) 



Due to the vertical direction of the magnetic field in arctic 
regions, the current density has only x and y components 

J = o {VyBl - v^j) . ( 92 ) 

Make a change of variables for the vertical direction 

z+h=»z; dz=*dz; range D^<z<D =» D^+h < z < D+h. ( 93 ) 



Then 



Jx R = oB[-v^z£ - V zj + {v y+Vjpc) 



( 94 ) 



The Biot-Savart law then gives 

dB' = ^ ° ah ~ ^ 



4 Jt 



4 ti 



dxdydz. (95) 



(x^ + + z^) ^ 



Integration over x,y,z gives the components of the total 
induced magnetic field. 



\l „ 0 B r z-D*h - Vq 






(x2+y2 + z2)3/2 



® ^ dxdydz. 



The y-component equation is similar. The z-component is 



JJJ 4n Jx.-J, — »r.D..A 



e 



“iwc 



dxdydz, (97 ) 



Performing the integrations [see Appendix B] yields 
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( 98 ) 



B' = 






- i<*) t 



2k^\Jkl*kx ^to o+o)' 



-icjt Q-^x^ [e ^>fio _ , 
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( 99 ) 



Bl 






2/ 






kji^. 



ky^jkl+k 



Approximating the geometry for an on-the-ice measurement 

f2 = 0; £>o = 0; D =»«>. ( 101 ) 






-2<j)C 



^x(y) 



2 ^ 0 ) 0 +^^ k^sjkl^k 



x(y) 



(102) 



Bz = 



ip. 



!V^ 



,+ 0)" 



k^sjkl^kl ky^Jkl+k 



( 103 ) 



Letting D approach infinity greatly simplifies the following 
calculations. This procedure is justified due to the expected 
magnitude of the parameters 

k^ = ; D-lOO-lOOOm; =►6'^'^ < 0 . 002 < 1 . ( 104 ) 
A 100m 



These equations still have a k-dependence. To get the 
components as functions of only o and t, integrate over the k- 
dependence 
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1 



( 105 ) 



Bi, = 






L 



ic =«• 



2 ^Cl) 



dk. 



-\ioOBv^e-^‘^'^ 1 2n[ 

2^0) 0+0)2 J^o -^in 



( 106 ) 



where is the smallest wavenumber (corresponding to the 

longest wavelength) present in the x-direction. A limit on 

wavelength is realistic for the ocean, and required to keep 

the integral finite. Recall that k^ keeps the total power 

finite as k approached zero, k^ thus represents the point in 

the power spectrum at which the characteristic k'^ behavior 

begins, k^ is thus a reasonable choice for k„- . Set k„. equal 

to k . Then 
0 



ln(y^+l) 



( 107 ) 



Similarly for the y and z components 

ln(y/2+l) _ ln(y/2+l) (3I3B) 

2 ^ 0 ) \/o) 0+0)2 kj, 

For the z-component, integration was done only over the k 
components which gave rise to the respective term via the 
Biot-Savart formula. 



38 



The total induced magnetic field is now 




The horizontal component of the induced magnetic field is 



Again, the power spectrum was obtained by Fourier 
transform of the autocorrelation of B', which is equivalent to 
squaring the amplitude here 



D. MEASURED MAGNETIC SPECTRA 

The Arctic Internal Wave Experiment (AIWEX) was conducted 
in 1985 on the arctic ice at 74 °N, 145°W. In April of that 
year, P. Czipott and W. Podney used two induction coil 
magnetometers to measure the North-South horizontal component 
of the Earth's magnetic field. Fig. 6 shows the results from 
three runs. The spectra flattened out at low frequencies 
because the electronic gain declines at low frequencies. The 
region in the 10'^ -10'^ Hz decade was relatively unaffected. 
The spectra here showed a characteristic 1/f^ dependence. 




( 110 ) 



(ln(v^+l))2 

2 (oVw“) ko 



( 111 ) 



(In (v/T+l) ) ^ 
2 kl 



( 112 ) 
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They also showed an approximately hundred-fold magnitude 
variation between the three runs. Thus, large magnitude 
changes can occur over a temporal scale on the order of days. 
Note that the bottom spectrum basically coincides with the 
dashed line spectrum which corresponds to measurements made on 
land at Poker Flat, Alaska. 
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Pow«r ap«ctra for magnacomactr output £roo Runs 100 
(c«nt«r)/ 103 (top)/ and 107 (bottom) at th* ic« camp/ 
and Run 124 (dashad curva) from Pokar Flat. Tha scala 
is callbratad by tha magnatomater gain at 0.1 Hz; 
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fraquancy. Error bar shows tha 90% confidanca 
intarral . 



Figure 6: AIWEX Power Spectra for magnetometjer (Czipott, 

Podney, 1985) 
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IV. RESULTS 



A. EVALUATION OF INTERNAL WAVE MAGNETIC FIELD MODEL 

Numerical evaluation of the formulas from the model used 
the following values/parameters 

= AtzxIO-’’ — ; 0 = 4.0—; B = 25,000 nT; (113) 

4^=3 .33x10-3 ;^=1,05 xI0-3-£^; (114) 

dz m* sec 

^max = 10^ = 6 . 28x1 (H5) 

Table II lists sample numerical values using the continuum w- 
and k-dependence method. Table III lists sample numerical 
values using the modal method. For the continuum method. Fig. 
7 shows the AIWEX model spectrum, the CEAREX model spectrum, 
and an extrapolated spectrum for the AIWEX magnetometer data. 
Extrapolation was made from the 10*^ - 10'^ Hz range of Run 
100. Fig. 8 shows the same spectra using the modal method for 
the AIWEX and CEAREX models. The vertical extent, D, of the 
internal wave field was 300 m (D'Asaro, Morehead, 1991) . 
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B'^2 (nT^2/Hz) 




Figure 7: AIWEX Model Spectrum, CEAREX Model Spectrum, and 
Extrapolated Spectrum for AIWEX Magnetometer Data (Continuum 
Method) 
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B>-2 (nT"^2/Hz) 




Figure 8: AIWEX Model Spectrum, CEAREX Model Spectrum, and 
Extrapolated Spectrum for AIWEX Magnetometer Data (Mode 
Structure Method) 
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Table II: Internal Wave Model Numerical Values 

(Continuum Method) 



frequency (Hz) B' AIWEX (nT^Hz) B' CEAREX fnT^Hz^ 




1.6x10'^ 



1.3x10^ 

1.3x10^ 

8.2x10^ 

1.7x10^ 



6.0x10^ 

6.0x10^ 

3.8x10'^ 

4.9x10^ 



Table III; Internal Wave Model Numerical Values 
(Mode Structure Method) 



Frequency fHz) B* AIWEX fnT^/Hz) B* CEAREX fnT^/Hz^ 



The AIWEX magnetometer measurements were not compensated 
for ionospherically generated magnetic fields. To separate 
the internal wave induced magnetic fields from the ionospheric 
fields, Czipott and Podney placed a tiltmeter on the ice as an 
independent measure of internal wave activity at the AIWEX 
site. They also obtained magnetograms of ionospheric activity 
during the experiment from several measurement stations in 
Alaska. Following Czipott and Podney, the signals from these 
three sources are compared in Table IV. 




2.7x10*^ 
2.7x10^ 
1. IxlO'^ 
1.2x10^ 



1.3x10® 

1.3x10® 

5.2x10^ 

3.3xl0’ 



1.6x10'^ 
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Table IV; Signal Characterizations 



RUN # 


IONOSPHERE 

(MAGNETOGRAM) 


TILTMETER 


MAGNETOMETER 


100 


"quiet" 


"high" 


"active" 


103 


"active" 


"active" 


"high" 


105 


"active" 


"quiet" 


"quiet" 


107 


"quiet" 


"quiet" 


"quiet" 



As a first order analysis, assume the ionosphere had only two 
levels of activity, "quiet" and "active", measured by the 
magnetograms. Assume the under ice seawater motion also had 
only two levels of activity, "high"/"active" and "quiet", 
measured by the tiltmeter. Then the above four runs show the 
resulting magnetometer signal for the four possible 
combinations of ionospheric activity and seawater motion. 
Runs 100 and 105 indicate that the ionosphere and the seawater 
motion each produced a measurable magnetometer signal. 
Together, runs 100, 103, and 105 indicate that the magnitudes 
of the magnetometer and ionospheric signals were of the same 
order. Run 107 affirms that the magnetometer readings result 
from actual ionospheric and seawater motion signals. 

Extrapolation of the AIWEX magnetometer data with a 1/f^ 
frequency dependence to the 10'^ - 10'^ Hz region gave power 
spectrum values between 10^ and 10^ nTVHz . For the continuum 
method, the model values varied one to two orders of magnitude 
above and below these extrapolated values. The l/k^jn 
dependence made these continuum models sensitive to the value 
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of chosen. Choosing = 1000 m instead of 10,000 m 
reduced the values by an order of magnitude. For the modal 
method retaining only the lowest order mode, the AIWEX model 
predicted values one to two orders of magnitude higher than 
the extrapolated data, while the CEAREX model predicted values 
about two orders of magnitude above the extrapolated data. 
The mode structure model values depended on vertical extent of 
the internal wave field as D^. Taking D = 100 m instead of 
300 m reduced the spectra by one order of magnitude and 
brought them within one order of magnitude of the extrapolated 
data. The vertical velocity profile of Equation [81] has a 
simple sin(5z) dependence. Refining this profile could 
significantly reduce the energy in each mode since is 
proportional to (velocity)^. Overall, the models reasonably 
predicted the magnitude of the seawater induced magnetic 
fields . 

The frequency dependencies of the continuum models varied 
significantly from the 1/f^ dependence of the extrapolated 
data. The AIWEX model predicted a linear increase with 
frequency and then a rapid falloff as the buoyancy frequency, 
N, is approached. The CEAREX model predicted a parabolic 
falloff with frequency as (N^-w^) . Unfortunately, the gain 
falloff of the AIWEX data over internal wave frequencies may 
have masked any dropoff near the buoyancy frequency. These 
continuum models predicted too much energy at the higher 
frequencies. Apparently the assumption of frequency and 
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wavenumber continue lead to erroneous predictions. The 
specific relation between o and k, denoted by in the 
models, must be determined. If one assumed that m(“) goes as 
1/ (number of wavenumbers), and that number of wavenumbers was 
proportional to frequency (see, for example, Garrett and Munk, 
1972, for a similarly modeled dependence), then the models 
would more closely predict the dependence of the data. One 
could compare the internal wave model in the 10'^ - 10'^ Hz 
range because internal wave oscillations are limited to 
frequencies below the buoyancy frequency. 

The mode structure method gave models with a much more 
realistic frequency dependence. The AIWEX model declined as 
1/f before going to zero near the buoyancy frequency. This 
dependence paralleled the observed AIWEX velocity spectra 
frequency dependence. The CEAREX model fell off as 1/f^, just 
as the CEAREX velocity spectra did. These parallel frequency 
dependences between induced magnetic field and generating 
velocities satisfied physical intuition. This close 
correspondence between model and data implied that a proper 
description of internal wave induced magnetic fields requires 
utilizing the internal wave modal structure. 

B. EVALUATION OF TURBULENCE MAGNETIC FIELD MODEL 

Numerical evaluation of the formulas use the following 
values/parameters 
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; B = 25, OOOnr; ( 116 ) 
(0 „ = 2tzx10"^Hz. ( 117 ) 



The value was determined the same way as for the 

internal wave model, assuming that the CEAREX data 1/f^ 
dependence extrapolates to the lO'^-lO'^ Hz range. Sample 
numerical values using Equation 102 are given in Table V. 
Fig. 9 shows a plot of the B„ ' power spectral density (nT^Hz) 
versus frequency for = 100 m and also the extrapolated 

AIWEX magnetometer spectrum. 

Table V: Turbulence Model Numerical Values 



FREQ (Hz) 


B'H^(nTVHZ) 


10'^ 


4.4x10''^ 


10'^ 


4.4X10'® 


5x10’2 


1.8x10'’ 


lO'’ 


4.4X10'’° 



(^n«x=100 m) 



The model spectra exhibited a characteristic l/o^ 
dependence, just as the AIWEX data show in the 10*^ - 10*^ Hz 
range. However, the model values were seven orders of 
magnitude less than the AIWEX data. At 10'^ Hz the model 
predicted magnitude was approaching the limits of the most 
modern field detectors at iC^ (nT^Hz) . The smallness of the 
predicted signal derived from the small velocities used and 
the maximum characteristic wavelength of 100 m. The small 
velocities on the order of 10’’ - 10'^ (cmVs^-Hz) were 
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experimentally measured during AIWEX (McPhee, 1989) and their 
power spectrum is shown in Fig. 10. The assumption of 
isotropic turbulence in the model restricted wavelength values 
to less than about 100 m because the maximum observed depth of 
the upper boundary layer is about 100 m (Stanton) (note: the 
"average" depth is about 30-40 m) . Thus, while turbulent 
motion in the upper boundary layer of the ocean could produce 
magnetic signals, the magnitude remained small because the 
volume of water generating the signal was limited. This 
implied that the AIWEX magnetometer measurements in the 10'^ - 
10'^ band were due to ionospheric signal only. 
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B’ (nT-'2/Hz) 



1.0E+09 
1.0E+08 
1.0E + 07 
1.0E + 06 
1.0E + 05 
1.0E + 04 
1.0E + 03 
1.0E + 02 
1.0E + 01 
1.0E + 00 
1.0E-01 
1.0E-02 
1.0E-03 
1.0E-04 
1 .OE-05 
1.0E-06 
1.0E-07 
1 .OE-08 
1.0E-09 
1.0E-10 




Turbulence 
AlWEX 



0.00001 0.0001 0.001 0.01 

Frequency (Hz) 



0.1 



Figure 9: Turbulence Model and Extrapolated 

Magnetometer Spectra 



AIWEX 
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log f s“ ^ 

-3 -2 -I 

T , ^ 1 . 1 T 1 ^ ^ ; 1 




2B.21 
29. 9B 

33.71 

33.^6 



Power sp>ecira of the venica] velocity time series at each level. The 
95 percent confidence interval is shown at lower left. The spectra are 
plotted on log-log axes, w’ith successive spectra displaced one decade down 
for clarity. 



Figure 10; AIWEX Vertical Velocity Power Spectra 
(McPhee, 1989) 
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V. 



CONCLUSIONS 



The magnetic fields induced by the motions of seawater in 
the earth's magnetic field were modeled analytically. Such 
analyses required mathematically modeling the velocities of 
the internal wave field and internal turbulence field, and 
then applying Ampere's law via the Biot-Savart formulation. 
For internal waves, the values predicted by the model for the 
induced magnetic field were roughly equal in order of 
magnitude to extrapolated magnetometer data. Treating 
internal wave frequency and wavenumber as continua yielded 
unphysical frequency dependences. Using an eigenmode 
structure for internal waves yielded more realistic frequency 
dependences. Thus, proper analysis of internal wave magnetic 
signatures must incorporate the internal wave modal structure. 
Additional magnetometer measurements should certify the 1/f^ 
frequency dependence and determine the behavior near the local 
buoyancy frequency. The turbulence model successfully 
coincided with magnetometer measurements in frequency 
dependence. The predicted frequency dependence matched the 
1/f^ dependence of the data. However, the predicted 
magnitudes were several orders of magnitude below the AIWEX 
magnetometer measurements and approached the limits of the 
most sensitive magnetometers. This implied that the 
relatively small scale of turbulence motions = 100 m) 
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resulted in a small induced magnetic field signature. The 
AIWEX data in the 10'’ - 10'^ Hz range were interpreted as 
deriving from ionospherically generated signals . 
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VI . RECOMMENDATIONS 



Additional measurements of seawater motion induced 
magnetic fields are needed. Such measurements should be made 
with highly sensitive magnetometers, for example, 
superconducting quantum interference devices (SQUIDs) . 
Discrimination against ionospherically induced magnetic fields 
must be done by using the large spatial coherence length of 
ionospheric fields to identify and remove their signal. 

The models can be extended to predict the signal seen by 
a gradiometer. This involves expressing the induced magnetic 
field at the field point in terms of field point position and 
taking the gradient. Gradient measurements of internal wave 
induced magnetic fields are more numerous than magnetometer 
measurements (for example, Czipott and Podney, 1985; Podney, 
1975) . 

Both the internal wave and the turbulence models for 
induced magnetic field utilized the approximately vertical 
orientation of the earth's magnetic field in the Arctic to 
simplify the calculations. The models can be extended to 
apply the Biot-Savart calculation method to the non-vertical 
magnetic field at other parts of the world's oceans. 

The internal wave model assumed a constant buoyancy 
frequency, N, for the entire wave field in the horizontal and 
vertical dimensions. In the actual ocean, the buoyancy 
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frequency changes from location to location and varies 
significantly with depth. A more refined model must include 
such variations and should predict a less defined cutoff 
frequency due to the spread in buoyancy frequencies. 

The modal version of the internal wave model should be 
refined by including the contributions of higher order modes. 
The resulting correction should not be too large because 
energy content rapidly declines for higher order modes due to 
dissipation. The velocity profiles should also be refined. 
As an example, require velocity continuity at the boundaries 
of the internal wave field and allow exponential decay outside 
these boundaries, analogous to going from an infinite square 
well solution to a finite square well solution. The amplitude 
within the well is reduced. 
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VII. APPENDICES 



A. CALCULATION OF INTERNAL WAVE INDUCED MAGNETIC FIELD 

For a region with constant Brunt-Vaisala frequency, N, the 
vertical component of the seawater velocity is given by 

w(x,y,z,C) = v(z) (118) 



where: 



6 






^ 2 - 0)2 



0 )‘ 



k*0 



( 120 ) 



Require as a boundary condition 

w{x, y, z=0 , t) =0 , => v(z=0)=l[q + C 2 ] , 



q=-C2 = C. 



( 121 ) 



Then the vertical velocity component becomes 

v(z) = . (122) 

The derivative with respect to z is 

-J-(v) = C(-iV^6) - C(-iic (123) 

oz ^ ^ 

With the simplified geometry, the y-component of velocity is 
zero. The x-component of velocity becomes 
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( 124 ) 



u(x,y, z, t) 



(i-^ - (4-) 

dz 



( 125 ) 



These expressions can now be used to calculate the induced 
magnetic field 



si . f f 

^ JxJz 2n r2 



( 126 ) 



+ « r Z 
Jx=-co 



z-D a (wB^ - uB^) 
271 i?2 



(h+z) dxdz 



( 127 ) 



(e“'-e-'“‘)co3Yco64> -f (e^'“*e-'“') 

r-rz-D B\x_^ 

J~m JzmD^ 2n (h+z)^ + 



(128) 

el'*-*-"') [h+z] dzdx 



The x-integration is done via contour integration in the upper 
half plane 




(h+z)^ + z^ 



dx 



{ 

Ju 






wfp (xti (^tz) ) (x-i (h+z) ) 



■ dx 



-L 






'bndry ( h + z) ^ + Z ^ 



-dx 



(129) 



i>Cjf(i (h+z) ) 

27ii[— 7— -] 
1 {h+z) +1 {h+z) 



+ 0 



71 -k^(h*z) (130) 

{h+z) 



Then 
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) cosYCOs<t) 



(131) 






i 






k. 



Performing the z-integrations and noting limits of integration 
yields 



2 



itit 



[ [( 












) e cosycos({) + 



032) 



+ r / 1 

-k^+ib 



^ (-kg*it) z_ 



k^-^ib 



7v) e 



(-Jc.-i6) z 



] ] Id„- (133) 



The z-component of the induced magnetic field is 
calculated similarly 

fx—cz-D \x^ di .. UB^) 



3' = /••*■*• r*'° iio^ X = r'i — 

* Jx.-» 2n J-» Jd, 2it 



X dzdx 



(134) 



= (e^‘--e-«^)cosYCOs<t. . asln(}> ^35 ^ 

J-Jd. 2n ' ^ (h*z)^ * 

= cosycos4>* 

2 II iCj, 

b' = ) e e cosycos<^ +(137) 

2 -ic^+i6 k*ib 






(138) 



Note that the x and z components have the same magnitude and 
differ in phase by w/2. These last equations can be 
simplified by multiplying the numerator and denominator by the 
complex conjugate, and using Euler's identity for complex 
exponentials, resulting in 

iBy=Bz- ^ — ) e'*''***’ e"-'"' [ (-^: sinSz - 6cos5z) 2icosYcos<t) + 

2 



+ (-k^cosbz + 6sin6z) 2&sin4) j 



(140) 
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B. CALCULATION OP TURBULENCE INDUCED MAGNETIC FIELD 

The integrations for the components of the turbulence 
induced magnetic field are shown below. 



/ H ^aB ry=<- rz=D*fi 



B' = 



^OJb rz-J^t 



(x 2 +y 2 + z 2 ) 3/2 



dxdydz, (141) 



\i^aB rx^<-rz^D*h -v^ /«y.. 2 dy dvdz 

TZ Jx—Jz-Do*h Jy*o (x^+y^ + Z^) ' ( 142 ) 



-2[i oBv rz^D*h ^x=- j . 

/ 2 I -7 r-r-j 7-^ dxdZ, (143) 



-2^ Qa5v^e“''“‘^ [■ z=r>^h z 

(x^ + z^) 



ATZ\jkl+klsj(ji\+(ji‘ 



^x=’>fz=D^h z e j j 

— 17 dxdz. (144) 

Jx=-«»^z=D„+?3 + 



Perform x-integration in upper half of complex plane 



B'x = 



4TZ)Jkl+k 



n dz, 

Jz^Dn-^h 



(145) 



V- o^Bv^ 



0---0 - e‘ 



2k^^kl+k 



. 



(146) 



The integration for the y-component is performed in the exact 
same manner, with x and y interchanged. 

The z-component integration is calculated similarly 



B 



,/ \i gCB ri<— rr— 

4 71 Jxm-"Jy^-mJz^D„*h 






y«o+fc>*^JCo+i:x (x^*y‘*z^)^^^ ^u>\*ut^^kl*ky (jf’+y^+z^) 



1(147) 
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4 7ly/<*) 



l-iMh ^ 1 ‘y.m 2xe^**^dxdydz . f^'“ 2ye^*^dxdydz ^ 

Wa Jx- Jy-O 7]^J7IJ(^2^y2.22)V2 V*Jx-0 ( (^2 ^ 2 ) 1/2 



2p„aBv„e-^"^ , 


r„».. 

z^D^*h J 


r x=*» 


X e^^^dxdz ^ , 


ry=- 


4n^ 




--V 




— y 



rz.p^A rx- x e^^^dxdz 



471^0^ 



' rz»D*n rx-- 
i/r-r.+h Jx--» 



/-. 



y e^‘^”^dydz 



yjkl*kl (x*iz) ix-iz) ^kl*ky (y*iz) (y-iz) 



2\x^oBv^e-^'^'^ rz-D.h 7iie'^>^ 



I* z-D+h ^ 

471^0)^0)2 \Jkl+kl \jkl*k 



] dz, 






2^0) 



+ 0)^ 



kjkl+k 



ky^kl+k 



( 148 ) 



(149) 



] (150) 



(151) 



(152) 
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C. INTERNAL WAVE VELOCITY POWER SPECTRAL DENSITY MEASUREMENTS 
Data from the Coordinated Eastern Arctic Experiment 
(CEAREX) , March-April 1989. 



Rotary Spectro 50m ADCP 




Frequency (cph) 
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Rotory Spectro 100m S4 




10 “' 10 ° 10 ' 10 -' 10 ° 10 ' 10 “' 10 ° 10 ' 



Frequency (cph) 
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qdo/^Cs /ujd) 



Rotory Spectra 1 50m S4 
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/s)Vcph (cm/s)2/cph 



Rotory Spectra 200m S4 



E 

<j 




«4 M« i . I < .i-.f ■« >*>»*!< — » « 4 »>>W) > i I 

10 "’ 10 ° 10 ’ 10 "' 10 ° 10 ' 10 "' 10 ° 10 ' 



F requency (cph) 
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(cm/s)^/cph (cm/s)^/cph 



Rotory Spectre 250m S4 
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D. 1-D TO 2-D POWER SPECTRUM CONVERSION CALCULATION 

The one-dimensional power spectrum and its derivative are 

1 






= -4-1. 



dk^ 






( 153 ) 



The two-dimensional power spectrum is given by 

dk^ 



n J k dk^ 

dk^ 



- k‘ 



__1 p 
" ki ^kl - k'‘ 



( 154 ) 

( 155 ) 



Make a change of variable 

/n = , K = -, dk^ = 

k^ ^ m ^ 



■ dm. 



ITT 



f^ik) = f 

^ 71 It 5 Jo 



4 p m*dm 

k^ Jo yiT 



m‘ 



( 156 ) 

( 157 ) 



From CRC 
M = a + 


Standard Math 
bm + cm^, then 


Tables 


(1987) , 


given a trinomial 


f —dm 


^ m"~^s/M - 


(2/3-1) h 


f dm 


- <"-l>^ f “'’■"dra. 




nc 


2nc 




nc J ^ 



Here , 
f,(k) 
f2(k) 



a=l, b=0, 



4 
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[ 



3 1 

4 k^ ■ 



c=-l. Applying this integral twice yields 

- ”Vl-'”" - * Isin-'m ]?, (159) 

4 8 8 ^ 

( 160 ) 
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